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THE  DISTRIBUTION  OF  THE  SIZE  AND 
NUMBER  OF  SHADOWS  CAST  ON  A  LINE 


The  present  paper  presents  algorithms  for  the  computation  of  the  distribu¬ 
tions  of  the  number  and  lengths  of  shadows  cast  on  a  line  segment  by  random 
obscurring  elements,  when  there  is  one  source  of  light  (in  the  origin),  and  the 
obscurring  elements  are  realized  according  to  a  Poisson  random  field*  These  al¬ 
gorithms  have  applications  in  problems  of  line-of-sight,  target  detection,  image 
processing  and  others. 

Key  Words:  Visibility  Probababilities,  Poisson  fields,  distributions  of  shadows. 

Acknowledgement:  The  present  research  was  partially  supported  by  the  U.S.  Army 
Research  Office,  Contract  DAAL03-89-K-0129  with  the  State  University  of  New  York  at 
Binghamton. 

1.  Introduction 

The  following  is  the  physical  problem  under  consideration.  One  source  of  light  is  located 
at  a  specified  point  0  in  space  (three  dimensional).  Consider  a  straight  line  C  of  distance 
r,0  <  r  <  oo,  from  0,  and  a  particular  segment  C  on  £.  Let  S  be  a  region  (strip)  in 
the  space  of  distance  u  from  0  and  distance  r  —  w  from  £,  where  0  <  u  <  w  <  r. 
A  countable  number  of  spheres  of  random  radius  are  randomly  dispersed  with  centers 
randomly  located  in  S.  The  location  coordinates  of  the  centers  of  the  spheres  and  their 
radii  constitute  a  Poisson  random  field.  These  random  spheres  cast  shadows  in  the  space. 
Some  of  these  shadows  fall  on  the  line  segment  C  and  obscure  it.  Due  to  the  randomness 
of  the  spheres  in  S ,  the  number  and  size  of  subiritervals  of  C  which  are  covered  by  shadows 
are  random  variables.  The  objective  of  the  present  paper  is  to  study  the  distributions  of 
these  lanuom  variables.  Chernorc  ' and  Daly  (1957)  studied  for  the  first  time  the  problem 
of  determining  the  distribution  of  the  length  of  a  random  shadow  on  a  line  segment  in 
a  Poisson  random  field  in  fhe  plane.  Explicit  algorithms  for  the  determination  of  this 
distribution  was  not  giver,  however,  in  that  paper.  In  a  series  of  papers,  Yadin  and 
Zacks  (1982,  1985,  1988)  developed  formulae  for  the  computation  of  the  probability  that 
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71  points  on  C  are  simultaneously  visible,  moments  of  random  weighted  visibility  measures 
and  approximations  to  the  distribution  of  the  number  of  shadows  on  a  line  segment  (19S4). 
Recently,  Yadin  and  Zacks  (1990)  applied  a  functional  which  was  introduced  by  Chernoff 
and  Daly  (1957)  to  derive  the  distribution  of  the  length  of  a  shadow.  In  the  present  paper 
we  extend  these  results  to  obtain  also  the  distribution  of  the  number  of  shadows  on  C. 

In  Section  2  we  introduce  basic  definitions  pertaining  to  the  random  fields  under  con¬ 
sideration  and  to  visibility  probabilities.  In  Section  3  we  present  the  distributions  of  the 
length  of  visible  and  of  shadowed  subintervals  of  C .  Section  4  is  devoted  to  the  recursive 
determination  of  the  probability  function  of  the  number  of  shadowed  subintervals  of  C. 
In  Section  5  we  present  numerical  examples  for  two  dimensional  and  three  dimensional 
models.  The  present  study,  as  well  as  the  above  referenced  ones  of  Yadin  and  Zacks,  were 
motivated  by  “lines-of-sight”  military  problems,  and  the  results  had  been  applied  by  oper¬ 
ational  analysts  dealing  with  simulated  war  games.  It  is  conceivable  that  the  methods  and 
results  of  these  studies  will  be  applicable  in  the  fields  of  image  processing,  when  portions 
of  the  image  of  an  object  under  consideration  arc  obscurred  by  randomly  dispersed  objects 
in  the  field  of  vision. 

2.  The  Random  Field  Model  And  Visibility  Probabilities. 

A  countable  number  of  spheres  (disks)  are  dispersed  so  that  their  centers  are  points 
in  a  region  (strip)  S  of  distance  u,0  <  u  <  oo,  from  the  origin,  0.  Let  C  denote  a 
random  vector  whose  components  are  the  coordinates  of  the  center  of  a  random  sphere 
(disk)  and  let  Y,  a  <  Y  <  6, 0  <  a  <  b  <  r  —  w,  denote  the  random  radius  of  such  a 
sphere.  Let  S*  =  S  x  [a, b]  and  B*  the  Borel  cr -field  on  S*.  Let  {(C;, Yi), i  =  1,2,...} 
represent  a  countable  sequence  of  random  disks.  The  random  vectors  (C,*,  Yi),t  =  1,2,. 
are  independent  and  identically  distributed,  having  a  c.d.f.  i7(c,y). 

Let  {B i, . . . ,  Bm } , 777  =  2, 3, ...  be  any  partition  of  S* ,  where  Bj  €  B*(j  =  1, . . .  ,m) 
and  let  N{Bj},  j  =  1, . . .  ,??r,  be  the  number  of  spheres  (disks)  having  centers  in  Bj.  The 
random  field  is  called  Poisson  if  N(Bj)  has  a  Poisson  distribution  with  mean 

(2.1)  fij  =  xj dH (c, y)  J  =  1,.. .  ,m,0  <  A  <  oo 

Bi 

ana  N~{B3 }  are  mutually  independent.  Let  C  be  a  segment  on  the  line  C.  A  point  on 
C  whose  distance  from  0  is  r ,  is  denoted  by  P0.  A  point  on  C  whose  distance  from 
Po  is  \x\  is  denoted  by  Pr,-^oo  <  x  <  oo.  x  is  the  coordinate  of  Pz  on  C.  A  point 
Px  £  C  is  called  visible  from  '0,  if  the  ray  Kx  from  0  through  Pz  is  not  intersected  by 
random  disks.  Let  A K±{x,i)  denote  the  expected  number  of  random  disks  centered  in  <S, 


2 


between  the  rays  R*  and  Rt  ,  which  do  not  intersect  1RX,0  <  t  <  co.  Explicit  formulae 
for  K±(x,  t)  will  be  gium  latci  for  some  special  cases.  Without  loss  of  generality,  assume 
that  the  intersection  of  5  with  the  plane,  £* ,  containing  0  and  £  is  a  strip  bounded  by 
two  parallel  lines,  U  and  >V,  of  distance  u,  and  iu  from  0;0  <  u  <  w  <  r  <  oo . 

Assume  that  the  segment  of  interest  on  £,C,  is  an  interval  whose  points  P*  have 
coordinates  x  €  [L,U].  Let  L*  <  L  and  U*  >  U  and  let  C*  be  a  region  in  S  whose 
intersection  with  the  plane  £*  is  the  trapez  bounded  by  Rj>  and  Ru* .  One  can  verify 
that  the  probability  that  Px  is  visible  is 

(2.2)  ip(x)  =  exp{-[/i{C*}  -  )'K-(x,x  -  £*)  -  -  x)}}. 

where  fi{C*}  is  the  expected  number  of  spheres  (disks)  having  centers  at  C* .  Formulae 
for  the  simultaneous  visibility  of  n  points  ,n  C  an  be  found  in  Yadin  and  Zacks  (1985, 
1988). 

3.  Distributions  Of  Length  of  Visible  And  Of  Shadowed  Segments- 
3.1  Distributions  Of  The  Length  Of  Visible  Segments. 

In  the  present  section  we  derive  a  formula  for  the  conditional  c.d.f.  of  the  length  of  a 
visible  segment  to  the  r.h.s.  of  Px,  given  that  Pr  is  visible. 

Let  I(x)  be  an  indicator  function  which  assumes  the  value  1  if  Px  is  visible,  and  the 
value  zero  otherwise. 

Let  L(x)  be  the  length  of  the  visible  setmenfc  of  C  to  the  r.h.s.  of  Px ,  i.e., 

(3.1)  L(x)  =  inf  {j/ :  y  >  x,  J(u)  =  1}  -  x. 

x<u<y 


We  derive  here  the  formula  for 

V(l  |  x)  =  P{L(x)  <  l\I(x)  =  1}. 

(3.2)  =  1  —  P'{l(x)  >  l\I(x)  =  1}. 

Let  C*  be  the  set  of  points  in  <S,  which  was  defined  in  the  previous  section.  We  derive 
the  ‘formuia  of  V^l\x) ,  for  L  <  x  <  U ,  and  0  <  T  <  u  —  x . 

Let  C~(x)  be  the  set  whose  intersection  with  the  plane  £*  is  bounded  by  Z/,  W  and  the 
rays  JR/,,  and  Rx.  Let  C(xJ)  be  the  set  whose  intersection  with  £*  is  bounded  by  U\  W 
and  the  rays  RX,JRX.}./;  and  (l  +  x)  the  set  whose  intersections  with  C*  is  bounded  by 
£/,W,lR/+x  and  Ry. .  Notice  that  C*  =  CT(.t)  U  C(x,  l)  U  C+(l  +  x) .  As  before,  we  denote 
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by  ii{C}  the  expected  number  of  disks  having  centers  at  the  sot  C.  Accordingly, 

P{L(x)  >  /,/(*)  =  1}  =  exp { — (//.{ C~(x ) }  -  \K-(x,x  -  L*)}  -  v{C(x,l)} 

-  [fi{C+(l  +  x)}  -  \I<+{1  +  x,  U*  -  l  -  *)]} 

(3.3)  =exp{-,/z{C'*}  +  A[tf_(a,*  -  L*)  +  K+(l  +  x,U*  -  /  -  »)]}. 

Dividing  (3.3)  by  (2.2)  we  obtain 

(3  4)  P{L(x)  >  l\I(x)  =  1}  =  exp{-A[24(x,  U*  -  x)  -  K+(l  +  x,U*  -l-  *)]}. 

3.2  The  Distribution  of  Shadow  Length. 

We  have  denoted  by  U(x)  the  right  hand  limit  of  the  shadow  on  C  to  the  r.h;s.  of  Px . 
Let  D(u\x)  denote  the  conditional  c.d.f.  of  U(x),  given  that  the  shadow  starts  at  Px. 

Consider  the  rays  Rx  and  Ry  for  y  >  x.  Let  N(xyy)  denote  the  number  of  disks 
centered  in  <S,  which  intersect  both  Rx  and  Ry.  Following  Chernoff  and  Daly  (1957),  we 
define  the  functional 

(3.5)  T(x)  =  sup {y  :  N(z,y)  >  1}. 

Furthermore,  let  Ti+1(e)  =  T(Tl(a;)),i  =  0, 1,  •  *  •  where  T°( x)  —  x.  Obviously,  T,+1(.r)  > 
T*( z),  for  all  i  >  0,  and  therefore  U(x)  =  Jim  T*(a;).  U(x)  —  x  is  the  length  of  the  shadow 

t— *  oo 

tothe  r.h.s.  of  Pz.  We  derive  first  the  c.d.f.  of  T(x).  Clearly,  {T(z)  >  t]  =  {N(xyt)  >  1}. 
Thus, 

(3.6)  P{T(X)  <  t}  =  P{AT(a;,/,)  =  0}  =  exp{-/i(rc,,i).}, 
where  ^(a;,<)  =  E{N(x,t)}  is  given  by 

(3.7)  (i{x,  t)  =  //{C7*}  -  A/4<®,  U*  -  a)  -  Atf_(i,  i  -  L*)  +  \K+feyt-  :r)  +  AJC.(*,  <  - 1), 

and  where  £  is  the  coordinate  of  the  bisector  between  Rx  and  Rt.  Notice  that,  since 
JCf(a;,0)  =  JC-(a;,0)  =  0  for  all  x , 

n{x,x)  =  lim//(a;,tf) 

U*  _ 

(3.8)  =.^{(?*}  -  A/4(s,-Cr  -  a)  -  \K-(xyx  ~  L*). 


Hence, 


which  is  the  probability  that  Px  is  visible.  Thus,  the  c.d.f.  of  T(m*i),H(t\x)  is  zero  for 
t  <  ft,  it  has  a  jump  point  at  x.Il{.r\x)  =  and  is  absolutely  continuous  at  t  >  x. 

This  property  is  inherited  by  the  c.d.f.  of  Tfl(a;),  Hn  (L\  x) .  We  provide  now  the  recursive 
relationship  between  Hn(t\  x)  and  Hn~i(t;x).  Introduce  the  bivariate  distribution 

On(t  -  P{Tn-\x)  <  tUT\x)  <  t2). 

Since  {Tn(x)  <t}c  {Tn~l(x)  <  *}, 

=  P{Tn(:.')  <  t) 

(3.10)  =Gn(t\t)  x\  all  t*>t. 


For  x  <  z  <  y  <  t) 

(3.11)  P{THx)  <  t\Tn~2(x)  =  z,Tn~l(x)  =  y}  =  exp {-[^y,t)  - 

Indeed,  given  that  {Tn~2(x)  =  z,Tn~l(x)  =  r/} ,  {Tn(x)  >  t]  if  and  only  if,  there  exists 
at  least  one  disk  which  intersects  Ry  and  R*,  but  does  not  intersect  Rr.  Hence, 

(3.12)  Gn(tl,t2\x)  =  J  J  exp{-(/i(w,<2)  -  n(z,t2)]}dG„-i(z,ti;x). 

These  bivariate  c.d.f.  can  be  determined  recursively,  starting  with  Gi(tiyt2\x)  =  H{t2\x) 
for  all  t\  <  t2.  Moreover, 

(3.13)  G7{h,t2-,x)  =  j'  c-',!-u'h)  Q*  e’^'^dz'j  dH(u\x). 

Finally,  since  Hn+i(t)  x)  <  Hn(t\x)  for  each  t  >  x  and  all  n  =  1,2,  •  •  •  the  c.d.f.  of  U(x) 
is 


(3.14) 


P{U(x)  <  t)  =  lim  Hn{t\x). 


Thus,  P{U(x)  <  t}  =  0  for  all  t  <  x,  and  lim P{U(x)  <  t}  =  The  conditional  c.d.f. 

tlx 

of  U(x),  given  {/(s)  =  0}  is 
(3.15)  P{V(x)<t\I(x).=  0}  = 


P{^x}<Pj  igs) 
1  —  ^(s) 

0, 


for  t  >  x 
for  t  <  x 


We  are  interested,  however,  in  the  conditional  c.d.f.  D(.u\?:) ,  where  Px  is  the  first  point 
(the  left  hand  limit)  of  the  random  segment  in  shadow. 


Simple  geometric  considerations  yield  that  the  length  ot  a  random  shadow  cast  by  a 
single  sphere,  having  left  hand  limit  at  Px,  which  intersects  C*  with  a  disk  with  center 
on  a  line  parallel  to  U  at  a  distance  h  from  O,  and  disk  radius  2,  is 


(3.16) 


U(xyh%z)  =  r  tan  2  sin 


■F  tan' 


-1  _ 

\rJJ 


where  (zC)  h )  are  the  coordinates  of  the  center  of  the  disk,  with 


(3.17) 


se=  Jfc  +  adi+0 


Thus,  if  a  <  Z  <  b  w.p.  1,  the  minimal  length  of  shadow  starting  at  x  is 

(3.18) 


Vm(x)  =  r  tan  (  2sin>»  \  (£)  (  !  +  (  *-  +  £  (l  +  g)’ 


+  tan"1  g).j 


Finally,  since  a  shadow  starting  at  Px  ends  at  point  U(x)  >  um(z)  •+ 

(3  191  D(u  I  x)  -  P{U ~  111  ~  P^U(-X)  ~  +  for  «  "  r 

(3.19)  D{u  |  x) - f_  P{U(7)  <  ^(X)  +  >  fol  u  *  x- 


4.  The  Distribution  of  The  Number  of  Shadows 

Let  J(x,y)>  for  L*  <  x  <  y  <  £/*,  denote  the  number  of  shadows  (invisible  segments) 
on  the  interval  (x,y)  on  C.  Let  Pj(x,y)  =  P{J(xiy)  =  j  |  I(x)  =  1}.  Obviously, 

(4.1)  P0(x,y)  =  1  -  V(y  -  x  |  x) 

The  stochastic  process  {I(X(t))yL*  <  t  <  U*}  is  a  regenerative  process.  A  cycle  in  this 
process  is  a  pair  of  visible  segment  followed  by  a  shadowed  segment.  Let  £(a)  denote  the 
length  of  a  cycle  starting  at  Pa .  The  distribution  of  ((a)  is  given  by  the  convolution 

px-a 

(4.2)  F(x\a)=  /  D{x  ~  a  |  t)dV(t  |  a). 

Jo 

It  fallows  immediately  that  the  probability  function  P3(x>y)  satisfies  the  recursive  equa¬ 
tion: 

(4.3)  Pj+1(x,y)=  f  Pj(t  +  x  |  y)dF{t  \  x),j  >  1 

Jo 


G 


and 


(4.4) 


Pik r‘>y)  =  /  [1  -  v(y  -t-x\t  +  x)]dF(t  I  x) 

Jo 

+  t'~X[l-D(y\t  +  x)]dV{t\x). 

Jo 


In  the  following  section  we  present  an  example  in  which  the  various  distributions  are 
approximated  numerically. 


5.  .An  Example 

In  the  present  example  we  restrict  attention  to  the  visibility  problem  in  the  plane.  The 
region  S  is  a  strip  with  boundaries  U  and  W,  parallel  to  £.  C  is  an  interval  on  £,  i.e., 
C  =  [:L,  U]  and  C *  =  [£*,  U*},  -oo  <  L*  <  L  <  U  <  U*  <  oo.  The  distances  of  U,W 
and  £  from  0  are  u,tu,r,  respectively,  0  <  u  <  w  <  r.  The  Poisson  field  is  a  standard- 
uniform  one,  i.e.,  centers  of  disk:',  «•*...  %C2,.  •  •  are  uniformly  and  independently  distributed 
in  <S,  and  the  radii  of  disks  are  independent  of  their  center  locations  and  have  an. identical 
uniform  distribution  on  (a,  b).  No  disk  in  S  can  intersect  0  or  £,  i.e.,  ( b  —  a)  <  u  and 
r  —  w  >  (b  —  a) . 


5.1  The  Functions  K±(x,t)  in  the  Standard-Uniform  Case,  With  Uniform  Dis-. 
tribution  of  ‘'adii  on  (a,  6). 

Let  K+(x,tyz)  denote  the  area  of  the  set  bounded  by  the  line  £+ ,  the  ray  R  x+t,t  >  0, 
and  the  lines  U  and  W;  £*  is  the-  line  parallel  to  Rx,  oil  its  r.h.s.,  of  distance  z  from 
it.  This  is  the  set  of  all  disk  centers  between  Rx  and  Rx-h,  of  radius  Z  =  z,  which 
do- not  intersect  Rx.  In  order  to  simplify  notation,  we  assume  that  w  =  r.  In.  actual 
computations  we  substitute  xw/r  and  tw/r  for  x  and  t  in  the  formulae  given  below.  Let 
d  =  (s2  +  w2)1/2 .  Simple  geometrical  considerations  yield: 

zd 


tf+0M>z)  =  -r{y  <uj 


w2  —  u2 
2w 


t-  2- 


(5.1) 


w  +  u 
2 


+  J{“-T<U’}  \^iw~zd) 


where  I{A)  is  the  indicator  set  function,  which  asumes  the  value  1  if  A  is  true,,  and  the 
value  0  otherwise. 

Notice  that  K+(x,  £  z)  depends  on  x  only  via  x2  .  Symmetry  implies  that  K-.(—x}  £  z)  = 
K+(x,t,z)  =  K+(— x,t, z)  for  all  -^oo  <  x  <  oo.  Hence,  K+(x,t)  =  K~(xyt)  and  wc 
delete  the  ±  subscript  of  K.  Finally,  K(x^i)  =  J ${K(x}t,Z)}  with  respect  to  the  uni¬ 
form' distribution  of  Z  over  (a, 6).  Let  x\  =  tu/d  and  X2  =  tw/d.  The  function  R{x>t) 
assumes  the- following  forms: 
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(i)  If  6  <  xi , 


(5.2) 


W2  -  u 2 
2w 


t- 


d 


u  -\-w 


(a -I-  *)  . 


(ii)  If  a  <  xi  <b  <X2 


w2  —  tr2  /  a?i  —  a 


2w 


6  -  a  u  +  w  6 


(53) 


2tw 


t  b  —  x\ 


+  ~  (  -  x\)  + 


d 2 


b  —  a  b  —  a 


3(6  -  a) 


(63  -  sj)  . 


(iii)  If  a  <  X;x  <  X2  <  b 


K(x,t)  = 


w2-v?  f  xi-a  d  1  ,  2-  2O 
“ST-  V  T-a"  ~  ~  a  ) 


(5.4) 

(iv)  If  xi  <  a  <  b  <  x 2 , 


+  2L  (<v  "  *?->  +  3(63#’ " x?))  • 


/e  Kv  r//  ,v  d  r  , v  ,  (P(a2 +  ab  +  b2) 

(5.5)  K(x,  t)=---(a  +  b)  +  - L. 


(v)  If  Xi  <  a  <  X2  <  b 

tw  X2  —  a 


(5.6)  K(x>t)  = 


2  b  —  a  2 (b  —  a) 


d  (^-a2)’+gJ„,/f  .^-a3). 


6<tu(&  —  a) 


(vi)  If  a2  <  a 


(5.7)  K(z}i)  =  0 

5.2  Algorithms  For  Discrete  Approximations. 

In  the  present  section  we  consider  discrete  approximations  to  the  functions  :Hn{t\  a;),  Gn 
(ti  j  ^2  >  )  n  ^  2. 

For  a  given  integer,  M,  partition  the  interval  (x,y)  to  M  subintervals.  Accordingly, 
let  8  =  (y  -  s)/M,to  =  £  and  ty  ==  U  ='0,1,*  •  •  ,M. 

For  i  =  0,  •  *  •  M,  let 

(5.8)  Hi(i)  =  JT('<i;<o)  =  expt{-— -/x(^0 » *« )}  • 
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We  compute  next  the  function  G2(h j)>i>  j  =  0, ■ ■ •  ,  M;  which  is  a  discrete  approximation 
to  (3.13).  For  i  =  0,  •  •  •  ,  M  and  j  =  ?.,••  •  ,  M , 

t 

(5.9)  G2(i,j)  =  ^exp {-(p(tk,tj)  -  K^.O))}  ■  [JTi(fc)  -  Hi(k  - 1)], 

t= o 

where  Arj(— 1)  =  0.  Notice  that  62(0,;)  =  #i(0)1  for  all  j  =  0, 1, •••  ,M;  and  for 

i 

(5.10)  G2(iJ)  =  -ffi(O)  +  exp{n(t0,tj)}  ■  ^exp{-^(ijt7i,)}[ffi(fc)  -Hi(k-  1)]. 

k~l 

Moreover, 

(5.11)  G2(i,j)  =  G2(j,j)  for  alii  >  j. 

Next,  compute  recursively,  for  every  n  >  3,  i  =  1,  •  •  •  ,  M,  j  =  i,  •  •  •  ,  jW 

t  i 

Gn(«»j)  “  K4».-,<;))}[<?n^l(fc>0  ~  <?n-l(fc  ~  1,0 

A-=0/=^- 

(5.12)  -  Gn-\(k,  l  -  I)  +  Gn-^k  -  1,/  -  1)], 
and  for  z  >  j 

(5.13)  Gn(i,j)  =  G„(i,i). 


For  z  =  0,Gn(0,  j)  =  C?n-i(0,  j)  =  tfi(0),  j  =  0,* • •  ,M.  After  computing  these  functions 
we  determine  Hn(j)  =  Gn(i,j),j  =  0,1,* ••  ,M.  jyn(;)  is  the  discrete  approximation  to 
the  c.d.f.  oil  Tn(x),  namely  Hn(t;x)\  i.e.,  Hn(tj\x)  «  ^n(i). 

In  Table  5.1  we  present  numerical  results  obtained  by  applying  this  algorithm  to  the 
following  special  case. 

5.3  Numerical  Example. 

We  compute  the  numerical  example  for  Table.  5.1  with  .the --following  geometrical  pa¬ 
rameters:  r  =  100[m],u  =  40[m],tu  =  60[m],a  =  l[m],6  =  2;5[m],a;  =  10[m],jD*  = 

-  +> 

— 100[?n],77*  =  100[m].  We  present  in.  the  tables  the  values  of  if„(j),n  =  1,2,3,;  = 
0,  •  •  •  ,20-,  when  8  =  l[m) . 
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Table  5.1.  Values  of  Hn(j )  for  two  values  of  A. 


A  =  0.02[l/??r2 

A  =  0.2[l/m2] 

j 

HiU) 

Mi) 

^3  O') 

SiU) 

Mj) 

HzU) 

0 

0.4914 

0.4914 

0.4914 

0.0000 

0.0000 

0.0000 

1 

0.5434 

0.5434 

0.5434 

0.0000 

0.0000 

0.0000 

2 

0.6009 

0.6008 

0.6008 

0.0000 

0.0000 

0.0000 

3 

0.6645 

0.6642 

0.6442 

0.0003 

0.0003 

0.0003 

4 

0.7342 

0.7336 

0.7336 

0.0021 

0.0020 

0.0020 

5 

0.8039 

0.8030 

0.8030 

0.0130 

0.0119 

0.0119 

6 

0.8664 

0.8652 

0.8652 

0.0576 

0.0522 

0.0522 

7 

0.9182 

0.9170 

0.9170 

0.1828 

0.1663 

0.1663 

8 

0.9569 

0.9559 

0.9559 

0.4164 

0.3843 

0.3842 

9 

0.9811 

0.9804 

0.9804 

0.6842 

0.6484 

0.6482 

10 

0.9934 

0.9930 

0.9930 

0.8763 

0.8507 

0.8505 

11 

0.9985 

0.9984 

0.9984 

0.9707 

0.9590 

0.9588 

12 

0.9999 

0.9999 

0.9999 

0.9981 

! 

0.9958 

0.9957 

13 

1.0000 

1.0000 

1.0000 

1.0000  ; 

1.0000 

1.0000 

14 

1.0000 

1.0000 

1.0000  1 

1.0000 

1.0000 

1.0000 

15 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

As  seen  in  Table  5.1,  the  convergence  of  Hn(j)  to  the  c.d.f.  of  U(x)  is  quite  rapid.  We 
have  therefore  approximated  the  c.d.f.  D(u  |  x)  by  the  sequence  D(j  |  i)  =  D(tj  |  t;),z  = 
0,1,*  ••  ,M,j  =  i,i  +  1,***  .  The  c.d.f.  of  ((a)^F(x  |  a)  was  computed  by  numerical 
convolutions  of  D(y  |  t )  with  V[t  |  x).  Finally  the  probability  functions  Pj(r,y)  were 
computed  recursively.  In  Tables  5.2  we  present  the  functions  Pj(x,y)  for  the  interval 
(—30,70)  for  various  values  of  A.  The  field  parameters  are  r  —  100, u  =  40,  w  =  60, a  = 
0,6  =  3,P*  =  —100,6/*  =  100  with  6=1.  We  see  in  Table  5.2  that  the  number  of 
shadows  on  (—30,70)  is  stochastically  increasing  with  A  up  to  A  =  0.012.  For  A  =  0.024 
the  number  of  shadows  starts  to  decrease  stochasticalfy.  We  have  here  a  delicate  interplay 
between  length  and  number  of  shadows. 
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Table  5.2.  Probability  Distributions  of  J(x,y), 
for  A  =  0.003, 0.006, 0.012, 0.024;  x  =  -30,  y  =  70. 


j 

A 

0.003 

0.006 

0.012 

0.024 

0 

0.0490 

0.0024 

0.0000 

0.0000 

1 

0.1841 

0.0270 

0.0010 

0.0016 

2 

0.2925 

0.1147 

0.0277 

0.1863 

3 

0.2598 

0.2414 

0.1600 

0.4639 

4 

0.1444 

0.2840 

0.3156 

0.2713 

5 

0.0535 

0.2020 

0.2891 

0.0670 

6 

0.0137 

0.0925 

0.1476 

0.0091 

7 

0.0025 

0.0287 

0.0471 

0.0008 

8 

0.0003 

0.0063 

0.0001 

0.0000 

9 

0.0000 

0.0010 

0.0016 

0.0000 

10 

0.0000 

0.0001 

0.0002 

0.0000 
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